Thermal rippling behavior of graphane 
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Thermal fluctuations of single layer hydrogenated graphene (graphane) are investigated using 
large scale atomistic simulations. By analyzing the mean square value of the height fluctuations 
{h^) and the height-height correlation function H{q) for different system sizes and temperatures 
we show that hydrogenated graphene is an un-rippled system in contrast to graphene. The height 
fluctuations are bounded, which is confirmed by a H{q) tending to a constant in the long wavelength 
limit instead of showing the characteristic scaling law q^~^ (rj 0.85) predicted by membrane theory. 
This unexpected behaviour persists up to temperatures of at least 900 K and is a consequence of 
the fact that in graphane the thermal energy can be accommodated by in-plane bending modes, i.e. 
modes involving C-C-C bond angles in the buckled carbon layer, instead of leading to significant 
out-of-plane fluctuations that occur in graphene. 
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Hydrogenated graphene (GE), called graphane (GA), 
is a quasi two-dimensional (2D) structure of carbon (C) 
atoms ordered in a buckled honey-comb lattice covalently 
bonded to hydrogen (H) atoms in an alternating, chair- 
like arrangemenlpi. Experimentally, it has been shown 
that GA can be obtained reversibly starting from a pure 
GE layer^ and since then it has become a material of high 
interest due to its potential applications in nanoelectron- 
ics^. As compared to GE, the chemisorption of the H 
atoms is accompanied by an important reconstruction 
of the chemical bonds and angles in the flat honeycomb 
lattice^. Each carbon atom acquires an H neighbor, in- 
volving a transition from sp'^ to sp^ hybridization, which 
turn the conjugated, graphitic C-C bonds into single C- 
C bonds changing locally the planar shape of graphene 
into an angstrom scale out-of-plane buckled shaped mem- 
brane^ as displayed schematically in Fig. [H 
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FIG. 1: (Color online) (a) Schematic view of a large sample 
of graphane. (b) Buckling between C atoms in the A (higher) 
and B (lower) sublattices at T = 300K. 



One expects that at non-zero temperature thermally 
excited ripples will distort the lattice. For GE these 
ripples can be described theoretically by the elastic- 
ity theory of continuum membranes^. Amongst others, 
this membrane theory predicts a suppression of the long 
wavelength out-of-plane fluctuations by the anharmonic 
coupling between the out-of-plane bending and in-plane 
stretching modes leading to a characteristic power law be- 
havior for the height fluctuations as a function of the sys- 
tem size. Although the height fluctuations still diverge, 
the normal-to- normal correlation is preserved over a large 
lenght scale stabilizing the membrane. By using atom- 
istic simulations these results were shown to be appHcable 
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to GE, implying that this prototype 2D solid behaves as 
a membrane^^^, and also to bi-layer GE^. Experiments 
have confirmed that suspended GE is not perfectly flat 
but instead presents ripples at finite temperatures^^. 

In this paper we study thermally excited ripples in GA 
using state-of-the-art molecular dynamics (MD) simula- 
tions and show that the angstrom scale thermal ripples 
present in GE do not appear in GA for temperatures up 
to at least 900 K. The A- and B-sublattice buckling is 
preserved and inhibits the formation of long wavelength 
ripples. As a consequence, the calculated height-height 
correlation function H{q) does not follow the character- 
istic q^~'^ power law scaling in the long wavelength limit 
predicted by membrane theory, and the height fluctua- 
tions appear to be bounded. 

According to membrane theory and within the har- 



monic approximation, applicable in the short wavelength 
regime where q is larger than some cross-over value g*, 
the out-of-plane and in-plane modes are decoupled and 
the elastic bending free energy density is described by 
Fharm = /^(V^/i)^, whcrc h is the local height and n is 
the bending rigidity of the membrane which governs the 
properties of the temperature induced ripples. Substitu- 
tion of the Fourier transform of h and integrating over 
2D space leads to the following height-height correlation 
function 



Hharm{q) = mq)^') 
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where N is the number of atoms, the surface area per 
atom and /c^ is the Boltzmann constant. Accordingly, 
the height fluctuations in the harmonic regime behave 
as {h^)harm = CL'^^ with C a temperature dependent 
constant and L the linear size of the system. In the large 
wavelength limit, i.e. for q < q^ ^ the height fluctuations 
are suppressed by the mentioned anharmonic coupling 
between bending and stretching modes giving rise to a 
renormalized g-dependent bending rigidity f<iR oc q~^ and 
a power law scaling behavior 
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and accordingly {h'^) = C'L^ ^ with C a constant 
(7^ C). The universal scaling exponent r] has been es- 
timated to be 0.821^^ For GE, usins; MC simulations 
with the empirical LCBOPII potential^^, good agreement 
with these results, derived from continuum theory, was 
found with 77 ^ 0.85 [7]. Here, we investigate to which 
extent membrane theory can be applied to the descrip- 
tion of thermally excited ripples in GA and we compare 
our results with those for GE. 

To calculate the height fluctuations for GA we flrst 
need to deflne an appropriate value for the height hi of 
each lattice site i. Since we are mainly interested in the 
long wavelength fluctuations, which normally govern the 
size of the height fluctuations, we deflned it on the basis 
of the carbon positions as 



hi = 



-E' 



where i is a carbon atom, Yl'j runs over the three carbon 
neighbors of i and Zi is the z-coordinate perpendicular 
to the plane. This deflnition allows for a straightfor- 
ward comparison with GE, for which the heights were 
deflned in the same way following previous work^. To 
measure {h'^) and H{q) we have performed MD simula- 
tions using the LAMMPS package^^. Both the GA and 
the GE systems have been sampled using the constant 
NPT ensemble (with P=0). For the interatomic interac- 
tions we used the modifled second generation of Brenner's 
bond-order potential, i.e. AIREBO^^, which has been 
shown to predict correctly the conflgurations for many 
different hydrocarbon structures. We simulated square 
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FIG. 2: (Color online) Heights of the C atoms in the GE 
(a) and GA (b) against the site index for arbitrary snapshots 
taken during the M P sim ulation. Ar=4860. T=300 K. (c) 
{h^) against L — y^LxLy in GE (circles) and GA (squares). 
H[q) for different system sizes as indicated for (d) GE and (e) 
GA. The dashed line shows the harmonic behavior and 
the solid line the correction due to anharmonic coupling for 
small q. Vertical arrows roughly indicate ^* below which the 
harmonic behavior is broken. 



shaped systems with the number of C atoms equal to 
N = 4860, 8640, 19440 and 37888, in a temperature range 
from 100 - 900 K (note that for GA the total number of 
atoms is twice large). Periodic boundary conditions were 
applied in the x- and y-directions. The presented results 
have been computed averaging over 300 — 400 uncorre- 
lated configurations. 

In Fig. [2] we show the out-of-plane positions for GE 
(a) and GA (b) for arbitrary snapshots taken during 
the simulation at room temperature. In GA, the A- 
and B-sublattices fiuctuate around their mean heights 
{hA,B) = ±0.2561. The scaling of {h'^) with the system 
size is displayed in Fig. ^c). The obtained values for 
GE are in close agreement with previous reported MD 
results obtained with the REBO potential^^ and slightly 
lower than those obtained from MC simulations using 
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FIG. 3: (Color online) (a) H(q) for different % of H atoms 
(0 % correspond to GE and 100 % to GA). The different 
curves have been shifted for a better comparison. In the inset 
we show the variation of (h^). T = 300i^. H(q) for different 
temperatures as indicated for (b) GE and (c) GA. The inset of 
(c) shows the average value of the heights in the A-sublattice 
of C atoms against temperature. N = 8640. 



the LCBOPlt^. For GE, (h^) increases as L^-^ as ex- 
pected from membrane theory. For GA, instead, {h'^) is 
almost independent of the system size. The differences 
in the rippling behavior of GA and GE are also evident 
from the results for H{q) (Figs. [2fd, e)). As it should 
be, the H{q) functions for different system sizes overlap. 
However, for GA, although the harmonic q~^ behavior 
for short wavelengths is well recovered, H{q) tends to a 
constant in the long wavelength limit. Hence, it does not 
follow the q^~^ power law as expected from membrane 
theory and found for GE. 

The intermediate regimes where GE is only partially 
covered by H atoms are analyzed in Fig.[3l^a). Notice that 
H{q) displays harmonic behavior even for a H covering 
as large as 90%. The deviations at low wave vectors are 
small and hence can be ascribed by anharmonic coupling. 
The variation of {h'^), shown in the inset, indicates that 
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FIG. 4: Variation of (a) (/i^), (b) q* and (c) k against tem- 
perature for GE (circles) and GA (squares). N — 8640. 



first the sheet is softened when partially hydrogenated 
and becomes stiff at full coverage. 

The behavior of H{q) for different temperatures is 
shown in Fig. [3l^b) and (c) . As expected, with increas- 
ing temperature H{q) is shifted to larger values for both 
GE and GA. In the inset of Fig. [3fc) we show the av- 
erage heights of the C atoms in the A-sublattice against 
temperature. We also found that {Ha) — —{hs) over 
the whole temperature range, implying that the A- and 
B-sublattice buckling is preserved. 

More signals of the differences between the corruga- 
tions in GA and GE comes from the temperature depen- 
dence of {h'^) shown in Fig. HJ^a). Note the difference in 
vertical scale displayed on the left- and right-hand side for 
GE and GA, respectively. While in GE the value of {h'^) 
changes about 0.61 A between 100 and 900 K the vari- 
ation is only 0.07 A in GA indicating that GA remains 
approximately un-rippled even at 900 K. The variation of 
q^ against T (Fig. [H^b)) also shows the same almost con- 
stant behavior for GA, whereas for GE, q*(T) behaves 
as expected for a 2D membrane^. From the calculated 
H{q oo) at different T and Eq. ([1]) one can also deter- 
mine the T dependence of the bending rigidity n. Using 
the REBO potential it was found that n decreases with 
T for GEi^, similarly as for liquid membranes. However, 
the opposite behavior was found using the LCBOPII po- 
tential^. In Ref. [l^ the rigidity was found to depend on 
the system size. Thus, the reported values for n in GE 
vary from 0.79 to 2.13 eV depending on the calculation 
method^^. In Fig. HJ^c) we show n for GE and GA calcu- 
lated from the harmonic part of H{q) between ^^=0.5 
and q=l confirming that for GE, n increased with 

T. For GA, n is much larger and most surprisingly, n 
strongly decays when temperature is increased, opposite 
to the behavior for GE. 

The reason why GA does not obey membrane the- 
ory should be found in the geometry of the buckled car- 
bon layer which allows for low energy in-plane bending 
modes, involving the C-C-C angles. These accordion- 
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FIG. 5: (a) Energy increase under bi- axial pressure with re- 
spect to the relaxed system at zero pressure for GA and GE 
(OK). Shape of (b) GE and (c) GA for P=9 MPa. Note that 
while GE is corrugated, GA instead remains un-rippled. 

type of modes have a relatively low energy. In principle 
this comes down to a strong anharmonic coupling of the 
out-of-plane bending mode with these in-plane bending 
accordion modes which strongly damp the out-of-plane 



excitations. To illustrate this further we have performed 
NPT simulations for increasing pressure, of which the re- 
sults are shown in Fig. [5l It shows that GA resists much 
higher pressures before bending than GE^. 

Using MD simulations we have shown that the intrin- 
sic thermal ripples present in GE do not appear in GA 
for temperatures up to at least 900 which we ascribe 
to the angstrom scale buckling of the carbon layer in GA 
into a carbon bilayer-like configuration. The rippling be- 
havior of GA is in disagreement with the continuum elas- 
ticity theory of membranes. The results from membrane 
theory are supposed to be universal which means that 
they should not depend on the atomic scale details within 
the membrane. Instead, we find that GA can accommo- 
date the thermal energy by in-plane bending modes, i.e. 
modes involving C-C-C bond angles in the buckled car- 
bon layer instead of leading to significant out-of-plane 
fluctuations that occur in graphene. The present results 
for GA suggests that membranes of atomic scale thick- 
ness can exhibit a more complicated behaviour than pre- 
dicted by membrane theory. 
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